function s = CoilCalibrationBellBloom(name)
    %% header
    import MeasurementSession.Magnetometer.MagneticFieldCalibration
    info = MeasurementSession.head;
    if nargin == 0
        info.session_name = 'MagneticFieldCalibration';
    else
        info.session_name = name;
    end
    
    p=mfilename('fullpath'); [~, idx] = find(p=='\');
    path = p(1:idx(end)); %path to this .m file
    
    %% Prepare
    s=MagneticFieldCalibration(info, [path 'device.dat'] );

    freqList=20:5:90;
    s.addParameter('BzScanFrom',      0, 'V');
    s.addParameter('BzScanTo',        8, 'V');
    s.addParameter('BzScanFreq',      0.1, 'Hz');
    s.addParameter('NIChannelIndex',  0:2);
    s.addParameter('NIChannelNames',  {'ScanVolt', 'LIA_x', 'LIA_y'});
    s.addParameter('NISamplingRate',  1e3, 'Sa/s');
    s.addParameter('AOMModFreqList',  freqList, 'kHz');

    %% Run
    s.run_measurement();

    %% Post treatment
    s.post_process();

    %% Finalize
    MeasurementSession.save(s);

end
